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ABSTRACT 



First results of a Cosmic Microwave Background (CMB) anisotropy 
experiment conducted at the Observatorio del Teide (Tenerife, Spain) are 
presented. The instrument is a four channel (3.1, 2.1, 1.3 and 1.1 mm) 3 Hc 
bolometer system coupled to a 45 cm diameter telescope. The resultant 
configuration is sensitive to structures on angular scales ~ 1° — 2°. We use the 
channels at the two highest frequencies for monitoring the atmosphere, and 
apply a simple method to subtract this contribution in channels 1 (3.1 mm) and 
2 (2.1 mm). The most intense structure at these two frequencies is the Galactic 
crossing with peak amplitudes of ~ 350 /iK. These crossings have been clearly 
detected with the amplitude and shape predicted. This demonstrates that our 
multifrequency observations allow an effective assessment and subtraction of 
the atmospheric contribution. In the section of data at high Galactic latitude 
we obtain sensitivities ~ 40 fj,K per beam. The statistical analyses show the 
presence of common signals between channels 1 and 2. Assuming a simple 
Gaussian auto-correlation model with a scale of coherence 6 C = 1.32° for the 
signal, a likelihood analysis of this section of data reveals the presence of 
fluctuations with intrinsic amplitude C^ 2 = 76^1 fiK (68 % CL including a 
~ 20% calibration uncertainty). Since residual atmospheric noise might still 
contaminate our results, we also give our result as an upper limit of 118 fiK at 
95% c.l. 



Subject headings: Cosmology: cosmic microwave background - Observations 



- 3- 



1. INTRODUCTION 

Measurements of fluctuations in the Cosmic Microwave Background (CMB) radiation 
provide one of the most direct tests of theories for the formation of structure in the universe. 
On large angular scales where the fluctuations are produced by the Sachs- Wolfe effect in the 
last scattering surface, the overall normalization Q rm s-ps ~ 18 fiK is well established by the 
COBE DMR flBennett et al. 19961 ) and the Tenerife beam-switching ( [Hancock et al. 1996|) 



experiments. On smaller scales (~ 1°) the acoustic effects on the last scattering surface are 
expected to enhance the level of fluctuations relative to those at larger angular scales, giving 
rise to a feature in the power spectrum of the CMB fluctuations, the so-called "Doppler 
peak" , with a position and shape which depends on the values of cosmological parameters 
(the curvature of the universe, the Hubble constant and the baryonic contribution). The 
current observations have recently started to allow a first determination of the position 
and height of this peak QNetterfield et al. 1996| , |3cott et al. 1996] ). Nevertheless more 



observations are needed to reduce uncertainties due to residual atmospheric signals or 
diffuse Galactic contamination. The main aim of the present experiment is to fill the gap in 
I space between the IAC-Jodrell Bank beam- switching experiments ( [Hancock et al. 1994 ) 



with a flat spectrum weighted I = 20 ±8, and both the ACME South Pole ( |Gundersen et al 
1995| ) {t = 68^32) and the Saskatoon experiment ( [Netterfleld et al. 1995| ). The experiment 



consists in a multichannel bolometer detector which is sensitive to multipoles in the range 
£ = 38 — 77 with the maximum sensitivity at I = 53. Here we present results of our first 
observing campaign in the summer of 1994. 



2. INSTRUMENTAL SET-UP 



The instrument is described in detail elsewhere ( [Piccirillo 1991| , |Piccirillo fc Calisse 
1993| ). In summary, the optics consist of a primary off-axis parabolic mirror (45 cm 
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diameter) coupled to a secondary off-axis hyperbolic mirror (28 cm diameter). The detector 
is a four channel photometer equipped with 3 He bolometers working at 0.33K. The bands 
are centered at 3.3, 2.1, 1.3 and 1.1 mm wavelengths (channels 1, 2, 3 and 4) as defined by 
a combination of resonant mesh filters. The instrumental noise is 3, 1, 1.6 and 1.2 mK- ^ sec 
in thermodynamic units for channels 1 to 4 respectively. High frequency leaks are blocked 
by a combination of fluorogold, black polyethylene and Pyrex glass filters. Laboratory tests 
have been done to check for the absence of significant leaks in the filters. The telescope 
is surrounded by 45° aluminum radiation shields fixed to the ground. The beam and 
side-lobes have been extensively checked by placing a distant Gunn source oscillating within 
Channel 1 band. An accurate bi-dimensional map of the beam shape has been obtained 
( 20° by 20°). We also studied the far field side lobe structures down to about —72 dB. 
These analyses show that the beam response can be approximated by a Gaussian with 
FWHM=2.°4 and that no significant side-lobes are found. 

The observing strategy consists of daily drift scans done at fixed position in azimuth and 
elevation. The beam throw in the sky is achieved by fixing the secondary mirror and chopping 
sinusoidally the primary mirror, according to 6{t) = 9 , (f)(t) = O + a x sm(2iif w t + ip) 
where 0(t) and 9(t) denote the azimuth and elevation respectively, (0 O , #o) = (0°,78.5°) is 
the initial position of the antenna (see next section), ao = 2.6° is the zero-to-peak azimuthal 
chopping amplitude at a reference frequency f w = 4 Hz. The demodulation of the signal 
is done online in software by evaluating the amplitude of the first (4 Hz) and second (8 
Hz) harmonic of the reference frequency. The resulting sky pattern for the transit of a 
point-like source resembles respectively the well known 2-beam and 3-beam response. In 
this paper we only deal with the second-harmonic demodulated data. We measured a stable 
DC offset (about 10 and 15 mK for, respectively, Chi and Ch2) which is partly due to the 
arc-shaped motion of the beam in the sky, i.e., the axis of rotation of the wobbling mirror 
is not exactly vertical, so the center position of the beam is a few arcminutes higher than 
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the two lateral position. The absolute pointing error has been estimated by observing the 
millimeter emission of the full Moon; conservatively we assumed that this error was not 
larger than the diameter of the Moon (~30'). The calibration constants, for each channel, 
have been chosen to give 1 K signal when the central beam is completely filled with a 1 K 
source. The system is calibrated using cryogenic cold loads and has been checked in the 
field performing hot/cold load tests, raster scans of the Moon and of the Gunn source. The 
calibrations are consistent to about 20 % absolute accuracy. 



3. OBSERVATIONS AND DATA PROCESSING 

The observations were carried out at Observatorio del Teide in Tenerife (Spain) 
during June and July 1994. This observing site is at an altitude of 2400 m and has been 
shown to have extremely good atmospheric transparency and stability (Watson et al. , in 
preparation). The precipitable water vapour during the campaign was below 1.5 mm for 10 
% of the time. The observed region of the sky was the strip at declination 5 = 40° where we 
collected about 550 hours of observations. This declination has been extensively measured 
from this site at larger angular scales (~ 5°) and lower frequencies (10, 15 and 33 GHz) 
with reported detections of structures in the CMB by the Jodrell Bank-IAC experiments 
( IHancock et al. 1994|) . 

The experiment was operated only during the night to avoid contamination from solar 
radiation. Fluctuations in the atmospheric emission are the main source of random noise 
in our system. The rms of the data collected during several consecutive hours shows that 
the sky noise of the data binned to 10 s was of the order of 5.1, 6.8, 11.5 and 13.6 mK 
for channels 1, 2, 3 and 4 respectively during typical observing nights. We experienced 
few nights with excellent observing conditions during which the noise dropped down by a 
factor ~3.5 in all channels. Even in these cases the atmospheric noise is larger than the 
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instrumental noise and therefore a major goal in our processing is to assess and subtract 
this unwanted source of noise. We can reduce the atmospheric noise in channel 1 and 2 
by subtracting the extrapolated signal from channel 4 (1.1 mm band) which is the most 
sensitive to atmospheric gradients. We checked also that extrapolating channel 3 (1.3 mm) 
produces similar results. In each channel i we have a superposition of astronomical signal 
and atmospheric signal: ATant,i = ^^'ant i + ^^ANTii where all terms are expressed in 
antenna temperature {Tant,i) and AT ant. i = ai AT ant a- Since the bulk of the signal in 
each channel is due to the atmospheric emission and the sky signal is expected to be much 
smaller, a linear fit of channel i versus channel 4 provides a very good estimation of on. For 
each channel we have to solve the following equation to obtain the sky signal: 



where fi and f$ are the atmospheric transparencies at channels i and 4 computed using 
measurements of the water vapour content, pressure and temperature of the atmosphere 
flCernicharo 1985j ). AT ant, i an d ATant,a are the data in channels i and 4 in antenna 
temperature units; AT^ stro is the astronomical signal in thermodynamic temperature units; 
Ci and C4 are the Rayleigh- Jeans to thermodynamic conversion factors:AT47VT.i = — ATi 
(cj =1.29, 1.66, 3.66 ,4.82 for channels 1 to 4 respectively) and pa is the fraction of the 
astronomical signals seen at channels i and 4: p i4 = 1 for CMB signal (i.e., outside the 
Galactic Plane crossing), pi^ ^ 1 in the area of the Galactic Plane crossing and evaluated 
according to our predictions for the Galactic foregrounds. Notice that the method requires 
a prior knowledge of the fraction of the signal expected in the channel i and in channel 4 
{pa). In the case of the crossing of the Galactic Plane, we use pa as obtained from our 
estimations. However, different models in the literature yield essentially the same values 
for the ratios pa and threfore we are confident that the recovered Galactic Plane is nearly 
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model independent. 

The atmospheric cleaning procedure was run on the data binned in 10 sec intervals, so 
that the noise is dominated by atmospheric noise in all channels. After cleaning, the 10 s 
binned scans are binned again to 4 minute bin size, so the beam is sampled with at least 
three points. The typical rms of these cleaned 4 minute binned scans are 0.33 and 0.29 mK 
for channels 1 and 2 respectively. These scans have an offset drift (always less than 0.6 mK 
for both channels) which we remove by fitting and subtracting combinations of sinusoidal 
functions with periods equal or larger than 72 degrees in RA, hence much larger than the 
scales at which our instrument is sensitive. We then proceed with obtaining the final data 
set by stacking all individual scans where the rms does not exceed 0.5 and 0.55 mK for 
channels 1 and 2 respectively. In this way we finally use 179 and 129 hours of data for 
channels 1 and 2 in order to form the final data set (33% and 23 % of the total data). The 
overlap in the data used in channel 1 and channel 2 amounts to 89 hours. Fig. 1 shows 
the results in thermodynamic temperatures obtained for channels 1 (a and b) and 2 (c and 
d) in the two sections of data away of the Galactic plane > 15°). The data have been 
binned in increments of l°in RA. The visual appearance of the results of channel 1 is nearly 
featureless, with all the points lying below the two-sigma level while channel 2 possibly 
shows the presence of signal. The mean error-bar in the results of each channel is ~ 70 fiK 
in a 1° bin in RA. In Fig. le and If we plot the sum ((Ch l+Ch2)/2) and in lg and lh 
the difference ((Ch 1— Ch2)/2) of both channels. 



4. GALACTIC FOREGROUNDS 



We have analyzed the Galactic contribution including the synchrotron, free-free and 
dust emission. The first two processes can be modeled by extrapolating the low frequency 
surveys at 408 MHz ( [Haslam et al. 1982| ) and 1420 MHz ([Reich fc Reich 1986|) . We use 
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a single power law to model the joint emission from these processes: Tff_ sync oc u^ 13 . In 
the Galactic Plane crossing at RA=305° the free-free emission is expected to dominate over 
synchroton due to the presence of the source Cyg A and other unresolved HII regions and 
therefore the spectral index is (3 ~ 2.1. This was additionally checked by obtaining the (3 
values from the 1420 MHz map, which reproduce the Galactic plane crossings at Dec=+40° 
both in the 33 GHz Tenerife beam- switching scans ( pavies et al. 19U6] , |Gutierrez 19920 
and in the 53 and 90 GHz COBE-DMR maps. We obtain /3 = 2.11 ± 0.05, which is the 
value that we adopted. With this value of j3 the free-free and synchrotron contamination is 
negligible at our frequencies. This is consistent with [Hancock et al. 1996| which found that 
at 33 GHz and 5 = 40° this contamination is less than 8 [iK on scales of about 5°. The 
dust contribution was estimated by using the 240 /zm DIRBE map as Galactic template 



and extrapolated to our frequencies using the model obtained in Poulanger et al. 1996 



According to this model the dust emission is well described by a grey body of emissivity 
7 = 2 ( [Reach et al. 19951) and with a dust temperature Td=17.5 K: l v oc u^B^Td). In 
Fig 2 we show the prediction (dashed line) for the Galactic emission and our results for 
channel 1 (top) and channel 2 (bottom). The general agreement between the predictions 
and our measurements constitute an important check on the performance of our system 
and our method of subtracting the atmosphere. Any reasonable combination of 7 and 
Td produces a negligible dust contribution outside the area of the main Galactic Plane 
crossing: rmsr,{chl) = 0.5 [iK and rmsn{ch 2) = 1.6 /j,K. In this paper we have ignored 
such contribution. A more detailed analysis will be presented in a forthcoming paper. 



5. STATISTICAL ANALYSIS 



We have analyzed statistically the data of channels 1 and 2 in the ranges 
RA 1= 241° - 285° and RA 2 =331° - 20° which are at Galactic latitudes \b\ > 15°. We 
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computed the correlation function to check for the presence of common structures in both 
channels. Fig. 3 shows the most relevant results obtained. In panel a the cross-correlation 
between channels 1 and 2 is shown, in panel b we show the auto-correlation of the combined 
scan (Ch 1+Ch 2)/2 and panel c shows the auto-correlation of the difference scan (Ch 
1— Ch2)/2. The error-bars in each of the three panels represent the 68% confidence 
levels and have been computed using Monte Carlo techniques. Due to the experimental 
configuration, the expected correlation in our data, from signals on angular scales at which 
the instrument is sensitive, will show a characteristic pattern with a positive feature at 
small angles, followed by a negative bump at angles ~ 4° — 10°, and then flat at larger 
angles. Visual evidence of a signal with this shape exists in both a and b panels, indicating 
a possible correlated signal between channels. The solid lines corresponds to the expected 
correlation for a model estimated using a likelihood analysis (see below). This model 
has an intrinsic amplitude of 86 /iK and a coherence angle 9 C = 2.1°, which corresponds 
to the amplitude and angle at which the maximum of the likelihood surface is attained 
for the combined analysis on channel 1 and 2. In c there is no evidence of structure, 
which is compatible with the expectations in the case of pure uncorrelated noise (dashed 
bands). Second, a likelihood analysis was used to determine the amplitude and origin of 
the signals. Our analysis assumed a Gaussian auto-correlation function (GACF) for the 
CMB signal: C intr {9) = C Q exp(-0 2 /20£). The GACF models together with its limitations 
and connections to more realistic scenarios have been widely discussed in the literature 
(see White fc Scott 1994Q and we adopt it to quote preliminary results, leaving physically 
motivated models for a future analysis. Table 1 presents a summary of the likelihood results 
for the coherence angle of highest sensitivity (9 C = 1.32°). Upper limits and detections are 
quoted at 95 % C.L. and 68 % C.L. respectively. The analysis of channel 1 does not show 

1 12 

evidence of signal in any of the two ranges considered with a limit C < 72 /iK, while for 
channel 2 we detect signal in both RA ranges. We made also a joint likelihood analysis 
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on channel 1 and 2 assuming that both channels have been measuring the same signal. 
This gives a detection C^ 2 = 76^21 A*K. Including the error in the absolute calibration we 
obtain C^ 2 = 76I32 /iK. Even though these results are compatible, we cannot exclude a 
possible frequency dependence of the detected signals. A two-dimensional (in the plane 

1 /2 

Cq-9 c ) joint analysis of both channels shows a well defined peak at C =86 /iK, 9 C = 2.1°; 
the expected correlation in our data for such a signal is plotted as the solid lines in Figs. 
3a and b. The likelihood analysis shows the presence of a clear signal in channel 2, while 
channel 1 is compatible with noise; this indicates a source of non-CMB contaminant in 
channel 2. It is unlikely that the data are contaminated significantly by the Galaxy (see 

1 /2 

previous section). If we convert our results from (C ,9 C ) to band power estimates, we 



find s J£(£ + l)C i /27r = 7.7tl\ 3 ■ 10~ 10 , marginally consistent with the Standard CDM 
model prediction of 1.7±8;| • 10~ 10 at our I = 53 QSteinhardt 19941) . The excess of signal 
seen in Ch2 might indicate that some residual atmospheric noise is still contaminating our 
results. Therefore we also quote an upper limit of 118 fj,K at 95% C.L. A more robust 
atmospheric subtraction technique together with a spectral index analysis of the measured 
fluctuations has been applied to our data set confirming that part of the signal seen in Ch2 
is of atmospheric origin. We will report these new results in a future paper. 
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RAi RA 2 RAi+RA 2 

Chi < 92 < 115 < 72 

Ch2 94tjg 1091^ lOllS 

Chi & Ch2 72152 79±il 76±t| 



Table 1: Table 1. Likelihood results in //K CMB. A 20% calibration uncertainty has been 
included. 
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FIGURE CAPTIONS 

Fig. I. — The stacked data sets at 3.3 mm (a and b), 2.1 mm (c and d), the addition (e 
and /) and the difference (g and h). In panel la we show the 3-beam profile indicating the 
instrumental response to a point source. Units refer to thermodynamic temperatures in this 
and following figures. 

Fig. 2. — Comparison of data at the Galactic Plane crossing at I ~ 80.2° (solid line) and 
predictions of foreground emission (dashed line) for channels 1 and 2 (see main text for 
details). 

Fig. 3. — Results of the correlation analysis: (a) the cross-correlation of Ch 1 with Ch 2, 
(6) the auto-correlation of (Chi + Ch2)/2 and (c) the auto-correlation of (Chi - Ch2)/2. 
The solid lines in (a) and (6) is the correlation obtained with a likelihood analysis, and the 
dashed line in (c) is the 95 % C.L. for uncorrelated noise. 
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